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Abstract 

Dynamics of eternal inflation on the landscape admits description in terms of 
the Martin-Siggia-Rose (MSR) effective field theory that is in one-to-one corre- 
spondence with vacuum dynamics equations. On those sectors of the landscape, 
where transport properties of the probability measure for eternal inflation are 
important, renormalization group fixed points of the MSR effective action de- 
termine late time behavior of the probability measure. I argue that these RG 
fixed points may be relevant for the solution of the gauge invariance problem 
for eternal inflation. 
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Consider a single field infiation driven by the infiaton field (p with potential V{(p). 
During one Hubble time At ~ the value of the infiaton field changes by 



A(f) ~ (j)At 



1 dV M|, dV 



3//2 d(b StiV 



under the action of the classical force — f^, and superhorizon fiuctuations 6(f) of the 
infiaton are generated with characteristic amplitude 



H 2V 

When 

M » A0, (3) 

the "quasiclassical" picture of the infiaton slowly rolling towards the minimum of its 
potential = (po, V{<Po) = breaks down, and we say that infiation enters eternal 
selfreproducing regime [1]. 

The physical picture emerging in this regime can be well understood in stochastic 
formalism [2,3]. As we know, the causal structure of eternally infiating Universe 
can be represented as an infinite set of Hubble patches with values of (p and the 
effective cosmological constant being random, stochastically distributed among 
the patches. The corresponding probability distribution P{(f),t) (its physical meaning 
is the probability density to measure a given value of the infiaton in a given Hubble 
patch) is the solution of the Fokker-Planck equation 

dmt) ^ ^ _J_d_ fdV \ 

dt 87r2 d(t)^ 3Hod(P\d(f) J ^ ' 

(for simplicity we consider V{(j)) = Vq + 6V{(f)), where |(^l^(0)| -C Vq; the Eq. (jl]) can 
be relatively easily generalized to the case of arbitrary potential V{(f}) [2]), and IR 
contributions to the correlation functions of the infiaton (including arbitrary loop 
corrections) can be simply calculated [4] by taking integrals 



m = J dcprp{M. (5) 

The general solution of the Eq. (jlj) is represented as the expansion 

P(0,t) = ^V^„(0)e-^"* (6) 



n=0 



over the eigenf unctions of the operator H. All eigenvalues En are non-negative since 
the "Hamiltonian" H is always supersymmetric. Moreover, if the normalization inte- 
gral 



' -iHi (7) 
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converges, the lowest eigenvalue of H is Eq = 0, so that at t — oo the distribution 
P{(p, t) approaches the time independent asymptotics given by 



-Pasymp(</') ~ exp 



3H^ J- 

This expression is known as the Hartle-Hawking wavefunction of the Universe [5] 111 

Discussing physics of eternal inflation, one was usually interested in this time- 
independent asymptotics because of two reasons: 1) time scale (5tasymp necessary for 
P to reach the asymptotics ([8]) is rather short in simple inflationary models: for 
example, for 6V{(f)) = im^^ we find (5tasymp ~ ^f^ = for 5V{(f)) = iA0^ — 
Stasymp ~ ^^^^^Hq^, etc, 2) the asymptotic distribution does not depend on 
initial conditions for eternal inflation, while the contributions of higher eigenstates 
into the expansion ([6]) do. 

However, in many physically interesting situations such as for eternal inflation 
on the string theory landscape [7, 8] the time scale ^tasymp is extremely long — very 
roughly one can estimate it as 5tasymp ~ ~ Nei^, where is the number of 
vacua on the landscape and ei is the eigenvalue corresponding to first excited state in 
the problem with a single dS vacuum. Also, often the time independent asymptotics 
(IH]) may not exist at all|^ In such situation, since it is clearly impossible to solve 
the Fokker-Planck equation (jl]) analytically for arbitrary V{(p), one can use various 
approximation schemes to analyze behavior of the probability distribution P{4>, t) at 
intermediate time scales Mp^ <^t <^ E^^. One such scheme is based on the method 
of dynamical renormalization group analysis [9], and it works especially well when 
transport properties of the landscape encoded in the Eq. (jlj) are important, such as 
for an infinite landscape with effective potential/Hawking- Moss/Coleman-de Luccia 
tunneling rates not growing at infinity [7, 8] . 

Very roughly, the idea is the following. Let us consider a n-point correlation 
function of the inflaton in the limit when all points are taken within a single Hubble 
patch. It can be represented in the following functional form 

m = J dci^rpiM = j d(t>r j -dps - hp^ (9) 

VPVPexp (i I dtd(f)P (dt - h) p) (10) 



since the probability distribution P{(p,t) satisfies the Fokker-Planck equation (jlj). In 
the limit t ^ oo (but t sufficiently small so that |P(0,t) — Fasymp(0)| ^ Pasymp{4>)) 



^It is interesting to note that the Linde-Vilenkin wavefunction of the Universe [6] is not a solution 
of the Eq. (jl]), so the stochastic formahsm allows us to unambiguously resolve the long debate 
between advocates of the Hartle-Hawking and Linde-Vilenkin wavefunctions. 

^This is the case for chaotic inflationary models with potentials V{(j)) = Xn4>", when the normal- 
ization integral ^ for the Hartle-Hawking wavefunction diverges at small </> [2]. 
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only renormalizable terms survive in the "effective action" 



5eff= / dtd(l>P{dt-H]P, (11) 



and scaling analysis of the latter allows us to learn many things about long time 
behavior of correlation functions (0"). In particular, it turns out that this behavior 
is determined by fixed points of the RG flow for the effective action fllip . Let me now 
present two examples of landscapes with RG fixed points. After that I will discuss 
how RG fixed points for the effective action ffTTj) may be relevant for the resolution 
of the famous gauge invariance problem for eternal inflation [10]. 

First, let us consider stochastic inflation in the potential with a weak disordered 
component y(0) =Vq + 51^(0), \5V\ < Vq [7]. In the limit \6V\ the solution of 
the Fokker-Planck equation (j4]) is given by 

F(0,t)~-^e (12) 

and we want to understand how the linear diffusion of the probability distribution 
P(0, t) is affected by weak disorder 5V{(f)) in the limit t ^ oo. 

For simplicity we take a gaussian distributed random potential with the correlation 
properties 

{5Vm = 0, (5^(0)5^(00) = A5(0 - 0') (13) 

where the averaging (■ ■ ■) is taken over the ensemble of potentials 5V , or, what is 
equivalent, over intervals (0, + 50) of one particular realization of 5V in this ensem- 
blei 

To determine how weak disorder affects dynamics of P(0,t), one can integrate 
it out in the "partition function" ([9]) and perform scaling analysis of the resulting 
effective action. One has 



-"eff 



IdtdS P^P-^P^P- 



(14) 



For qualitative understanding of physics it is actually sufficient to use simple scaling 
arguments: from the quadratic part of the effective action iHM . by multiplying 
by a factor of / we immediately find that the time scale t scales as P, the product 



■^Naively, it seems that the gaussian ensemble (|13p is not weh defined, because it contains con- 
figurations corresponding to negative values of the full potential V{(l)). However, if one changes the 
correlation properties (fT3|) accordingly to resolve this issue, one finds [8] that the change does not 
influence the long time behavior of the probability distribution P(0, t) since it introduces irrelevant 
terms into effective action (fTTj) . 
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PP — as and the disorder strength A — as as well. Therefore, we conclude 
that weak disorder of the random potential type fll3l) is irrelevant for the long time 
behavior of P{(f), t) and does not influence the linear diffusion law (0^) ~ t. The proper 
renormalization group analysis of the theory f|T^ confirms this conclusion [7,8]. 

Since we are also interested to know the intermediate time scale behavior of the 
probability distribution P{(j), t), we can try to figure out how important are the excited 
states ipk{4')^ k > present in the expansion for the dynamics of P{(f),t). As it 
turns out [8,11], for any given realization of the disorder all eigenstates are localized 
or peaked near the points (j) = 4>i corresponding to the deepest minima of the full 
potential V^((/>). The width of ipk{(p) near a localization center 0j behaves as 

growing with decreasing k. Excited, localized, eigenstates dominate in the expansion 
(E]) at small t, and one has 

{{C^ - CP.f) t\ (16) 

so the diffusion is slowed down significantly. At later times, the regime of diffusion 
becomes linear, as we have found above. 

Let us now turn to a somewhat different case of stochastic inflation in the presence 
of random force — when the correlation properties of disorder are 

'Z)-0, (17) 



After integrating disorder out we find the following effective action 

S,, = /*.m(p|p-^P^P- (18) 

Simple scaling analysis leads now to the conclusion that a scales as /, so that the 
disorder f[T7l) is important for the long time dynamics of the probability distribution 
P(0, t). More accurate renormalization group analysis shows that the RG fixed point 
cr = 0"^= 7^ exists [7,8, 12]. The latter determines dynamics of P(0, f) at t — oo, and 
the diffusion process is dramatically slowed down — instead of the linear diffusion 
law ((/)^) ~ t one has weak logarithmic behavior 

(02) ^ log^t. (20) 

Similar effect is present on the landscape of dS vacua with disorder [8], where "old" 
inflation becomes eternal due to the possibility of multiple consequent Coleman-de 
Luccia tunnelings between different de Sitter vacua. Dynamics of the probability 
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Pi{t) for an observer to find lierself in a bubble of dS vacuum with the cosmo logical 
constant Aj is described by the vacuum dynamics equations [13] 

dP 

^ = E (r.^^:;- - r,.,^.) = E H^^P^' (21) 

j 3 

where are Coleman-de Luccia tunneling rates (inverse characteristic times of tran- 
sition) between de Sitter vacua i and j. Technically, equations fl2T|) represent the 
discretized version of the multifield Fokker-Planck equation (jlj). Not surprisingly, it 
turns out that the case of weakly disordered tunneling rates [8] 

r,, = f % + 5T,j (22) 

is reduced to the two types of disorder (random potential and random force) discussed 
above. The effects of disorder are also similar: when the number of vacua adjacent 
to any given one is equal to 2, the diffusion of the probability distribution Pi{t) is 
suppressed as 

{n\t))r^\ogH (23) 

at t — ^ oo (here n is the index numerating dS vacua on the landscape). 

Let me now discuss how RG analysis of the effective action (ITT]) may help to resolve 
the gauge invariance problem for eternal inflation. The probability distribution P(0, t) 
governed by the Eq. (jl]) does not take exponentially inflating 3-dimensional volumes 
of Hubble patches into account. The reason for that is the effective coarse graining [2] 
at scale I ~ performed during the derivation of the Fokker-Planck equation (j4]). 
Usually, the effect of volume growth is taken into account [14] (see also [15] for recent 
developments) by adding the volume weighting term to the r.h.s. of the Fokker-Planck 
equation: 

dtPv = {^i[§jPy) + {H''\mH^'\ct^)Pv) + m<t^)Pv m 

Since H is positively defined, patches with larger 3-dimensional volume and higher 
H are clearly dynamically favored, which leads to many technical and conceptual 
problems associated with the measure Py(0, t). Probably, the most widely discussed 
one among them is the gauge invariance problem [10] — the probability Py and 
associated correlation functions of physical observables are exponentially sensitive to 
the choice of time reparametrization □ 

dt = f{T)dT. (25) 

Another important problem, directly related to the previous one is very strong depen- 
dence of the probability distribution Py on initial conditions for eternal inflation, as 

^One may even require stronger invariance of Py w.r.t. reparametrizations dt = /(r, (j))dT. 
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can be transparently seen for probability measures with effective volume cutoff [16]. 
Before I proceed to the discussion of RG analysis and RG fixed points of the effective 
action (1111) . a couple of remarks has to be made. 

First of all, one may observe that there exists a problem with Eq. even more 
serious than the problem of gauge invariance — the probability flow described by the 
modified Fokker-Planck Eq. fl2^ is non-unitary, i.e., the total probability J d(f)P{(p, t) 
is not conserved. It was understood long time ago [3] how write the Fokker-Planck 
equation for the normalizable probability — a non-local term proportional to {H) has 
to be added to the r.h.s. of the Eq. fl2^ : 

mM = a,(^^p(0,t)) + ^a,(//3/2(0)9,i/=^/2(0)p(0,t))+ (26) 

+3P(0,t)i/(0)-3P((/),t) j d<P'H{<P')P{(l)') (27) 

Dynamics of P described by the Eq. (!26l) is much healthier than the dynamics of Py- 
In particular, the exponential sensitivity of the probability measure P to the time 
reparametrizations (!25|) is lost, although the issue of gauge invariance is of course not 
yet resolved completely. 

Second, one actually has to ask herself — is it true that correlation functions 
of any physical observable should be invariant w.r.t. the reparametrizations fl25l) ? 
The answer to this question is certainly negative, and to illustrate this, the following 
counter-example can be presented. Instead of calculating correlation functions of 
the infiaton we can calculate correlation functions of the world time t = t{(f)) or the 
number of efoldings A^(0) = log a = J* Hdt'. The world time and the number of 
efoldings are two different quantities, related to physically different observables: t is 
essentially the quantum phase of a heavy (m ^ H) particle living in a given Hubble 
patch, while iV(0) is the phase of a light (m -C H) particle. Clearly, the correlation 
functions of these two observables should behave differently, and both observables are 
not invariant w.r.t. the transformations fl25l) . 

With these remarks in mind, let us still try to find whether solutions of the Eq. 
f l26|) exist which weakly depend on initial conditions and are invariant w.r.t. the 
transformations fl2Sl) . First of all, we note that volume weighting term 3HP and the 
non-local term —3{H)P cancel each other if time is measured in units of local Hubble 
parameter H, and one has 

^^"^■^^ = ^'^ {jHWf) ^ 8^^' {H'^\mH'^'i<P)P) ■ (28) 

The general solution of this equation has the form similar to ([6]). In particular, at 
log a ^ oo it approaches the Hartle-Hawking asymptotics ([8]). Of course, the latter 
a) does not depend on initial conditions since it is the attractor for the ensemble 
of trajectories P = P{(p,t) and b) is invariant w.r.t. reparametrizations of t simply 
because the expression ([H]) is the function of the infiaton potential V{(f)) only. 
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However, as was explained above, in at least two important classes of inflationary 
models this asymptotics is never reached because the Hartle-Hawking wavefunction is 
not normalizable: (a) infinite landscape of vacua a la condensed matter with effective 
potential/Hawking-Moss/Coleman-de Luccia tunneling rates not growing at infinity, 
when transport properties of the probability flow P become of especial importance, 
and (b) models of chaotic inflation where the regime of deterministic slow roll acts 
as a sink for the probability flow (jl]) and inflation ends after slow roll conditions 
break down. Fortunately, analogues of the Hartle-Hawking asymptotics exist for both 
classes of models which have the properties of invariance w.r.t. reparametrizations of 
t and independence of initial conditions. 

Let me start with the class of models (a), where the hint is provided by the 
dynamical renormalization group analysis of the effective action fllip . As we have 
discussed, long time dynamics of the probability distribution P{t) is determined by 
the trivial RG flxed point corresponding to the regime of linear diffusion (0^) ~ t 
unless non-trivial RG flxed points for the effective action ( fTT]) exist. These RG flxed 
points correspond to attractors for the ensemble of trajectories P = P{t). Once the 
probability distribution reaches such attractor, memory about initial conditions is 
lost@ On the other hand, P{t) is invariant w.r.t. time reparametrizations (!25|) due 
to the conformal invariance of the correlation functions at the RG flxed point. I have 
presented above one such example of the inflationary theory with dynamical RG flxed 
point: inflation on the low-dimensional landscape with disorder, when the behavior 
of the correlation function (0^) is logarithmic!^ 

Non-trivial stationary solutions of the Fokker-Planck equation also exist for chaotic 
inflationary models (class (b) above) [2]. To flnd them, one needs to rewrite the 
Fokker-Planck equation fl25]) in the form 



aioga 50' ' SnMf, ' Sir V 

Instead of the Hartle-Hawking wavefunction ([8]) the time independent asymptotics of 
the general solution of the Eq. (129!) is given by the solution with j = Const. Such 



^This is not so for the trivial RG fixed point corresponding to the hnear diffusion law (0^) ~ t 
— the gaussian distribution is always peaked at the same point where the initial distribution 
is peaked. 

®I would also like to mention another well known example of the theory with the dynamical RG 
fixed point, which perfectly illustrates why the analysis of RG fixed points is important — the theory 
of Kolmogorov's turbulence. The MSR effective action for the Navier-Stokes (or Euler) equation in 
the presence of stochastic noise has always at least two RG fixed points. One of them corresponds 
to the ultimate equilibrium, v = and is trivial, while another one is not — it corresponds to the 
Kolmogorov's cascade (of decaying turbulence). The latter is the stationary flow of energy from 
larger scales towards smaller scales. In the regime of cascade, correlation functions of observables 
(such as velocity of the fluid v) are scale invariant with non-trivial scaling properties, and any 
information about initial conditions for these correlation functions is lost as long as the regime of 
cascade is reached, since Kolmogorov cascade is the attractor. 
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solutions describe constant flow of probability towards the regime where inflation 
becomes deterministic (i.e., when the condition breaks down) and ends afterwards. 
They also have properties similar to the Hartle-Hawking asymptotics ([S]) — they are 1) 
dynamical attractors and therefore do not carry information about initial conditions 
and 2) transparently invariant w.r.t. the time reparametrizations fl25|) . 
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